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^ ■ We discuss the physics underlying event-by-event fluctuations in relativistic heavy ion colli- 
sions. We will emphasize how the fluctuations of particle ratios can be utilized to explore the 
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properties of the matter created in these collisions. In particular, we will argue that the fluctu- 
tions of the ratio of positively over negatively charged particles may serve as a unique signature 
for the Quark Gluon Plasma. 
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3 i 1. INTRODUCTION 

Any physical quantity measured in experiment is subject to fluctuations. In general, these 
fluctuations depend on the properties of the system under study (in the case at hand, on the 
properties of a fireball created in a heavy ion collision) and may contain important information 
about the system. 

The original motivation for event-by-event (E-by-E ) studies in ultra relativistic heavy ion 
collisions has been to find indications for distinct event classes. In particular it was hoped 
that one would find events which would carry the signature of the Quark Gluon Plasma. First 
pioneering experiments in this direction have been carried out by the NA49 collaboration |jlp. 
They analysed the E-by-E fluctuations of the mean transverse momentum as well as the kaon 
to pion ratio (see Fig. ([I])). Both observables, however, did not show any indication for two or 
more distinct event classes. Moreover, the observed fluctuations in both cases were consistent 
with pure statistical fluctuations. 

On the theoretical side, the subject of E-by-E fluctuations has recently gained considerable 
interest. Several methods to distinguish between statistical and dynamic fluctuations have been 
devised [0,0]. Furthermore the influence of hadronic resonances and possible phase transitions 
has been investigated [Q-|8|]. All these theoretical considerations assume that the observed fluc- 
tuations will be Gaussian and thus the physics information will be in the width of the Gaussian. 
In [|9|] one of us has shown, that in this case, the information contained in E-by-E fluctuations is 
that of a two particle density, which can alternatively be measured using a two arm spectrome- 
ter. This observation implies that possibly interesting fluctuations can be observed and verified 
by a whole array of detectors not only large acceptance ones. 
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Figure 1 . Results from the NA49 collaboration for the E-by-E fluctuations of the mean trans- 
verse momentum as well as the kaon to pion ratio [|T|] 



2. Event-by-event fluctuations 

Fluctuations have contributions of different nature. First there are the 'trivial' statistical fluc- 
tuations due to a finite number of particles used to define a particular observable in a given 
event. Second there are dynamical fluctuations. Those carry the information about the proper- 
ties of the system. In the special case of heavy ion collisions, there is another source of 'trivial' 
fluctuations, namely the fluctuations of the volume. All these fluctuations contribute to E-by-E 
fluctuations. The obvious challenge, therefore, is to separate the 'trivial' from the interesting 
components. 

2.1. Fluctuations in thermal systems 

In a thermal system, the fluctuations of a quantity are proportional to its susceptibility, which 
is the second derivative of the appropriate free energy with respect to the conjugate variable. 
For example the fluctuations of the charge are given by 

m = -Tp^- = VT XQ (1) 

where /jq is the charge chemical potential, T the temperature and V the volume of the system. 
Xq is the charge susceptibility. Here we will mostly concentrate on the charge fluctuations. 
Therefore, let us discuss the charge susceptibility in some more detail. First, as shown in []Hj], 
the charge susceptibility is directly related to the electric screening mass in a system. Further- 
more, it is given by the zero momentum limit of the static current-current correlation function 

m 

x e = n 00 (w = o,^^o) (2) 

which can be, and has been, evaluated in Lattice QCD JTT|JT2| ] as well as in effective hadronic 



models [10,13,14] 
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2.2. Fluctuations in Heavy Ion Collisions 

The systems created in a heavy ion collisions carry another important, but rather uninteresting 
source of fluctuations namely the fluctuations of the volume. Even for the tightest centrality 
selection the impact parameter and thus the volume of the system still fluctuates considerably. 
These volume fluctuations may even dominate. For example let us consider the fluctuation of 
the number of particles of a certain species. In a given event the number of particles is given by 

N = pV (3) 

where p is the particle density and V the volume of the system created in this event. Thus the 
fluctuations have contributions from both the density fluctuations and the fluctuations of the 
volume 

<(5^) 2 > = <(5p) 2 )<V 2 > + <p 2 ><(5y) 2 >. (4) 

The interesting physics of course resides in the fluctuations of the density. Consequently, one 
needs to consider so called 'intensive' quantities, which do no explicitely depend on the volume. 
One such observable is the ratio of particle abundances, and this is what we will concentrate on 
here. 

2.3. Fluctuations of particle ratios 

As discussed in the previous section, the fluctuations of particle ratios should be independent 
of volume fluctuations. This is certainly true if one looks at similar particles such as 7t + and 
where the freeze out volumes are expected to be the same. Some residual volume fluctuations 
may be present if one considers ratios of particles with different quantum numbers such as the 
K/n ratio, but they still should be small. Let us define the particle ratio R12 of two particle 
species N\ and N2 

The fluctuations of this ratio are then given by 

(ofli 2 ) 2 = f (m) 2 ) (m) 2 ) mm \ 
(Rn) 2 \ (m) 2 <Af 2 > 2 wmJ' 

The last term in Eq. (j^) takes into account correlations between the particles of type 1 and type 
2. This term will be important if both particle types originate from the decay of one and the 
same resonance. For example, in case of the 7t + /7t~ratio, the po, CO etc. contribute to these 
correlations. Also this term is responsible to cancel out all volume fluctuations [Q]. 

Let us note that the effect of the correlations introduced by the resonances should be most 
visible when (Ni) ~ (N2). O n me other hand, when (N2) ^> (Ni), as in the K to K ratio, the 
fluctuation is dominated by the less abundant particle type and the resonances feeding into it. 
The correlations are then very hard to extract. In [|5J] it was shown that in case of the K/n -ratio 
resonances and quantum statistics give rise to deviations from the statistical value of at most 
2 %, in agreement with experiment [[!]]. 

As pointed out in [|5|] the measurement of particle ratio fluctuations can provide important 
information about the abundance of resonances at chemical freeze out, and thus provides a 
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crucial test for the picture emerging from the systematics of single particle yields [JT5|] . In 
particular the fluctuations of the ft + /ft~-ratio should be reduced by about 30 % as compared to 
pure statistics due to the presence of hadronic resonances with decay channels into a 7t + -7t~-pair 
at chemical freeze out. About 50 % of the correlations originate from the decay of the po and 
the oo mesons. Thus the fluctuations provide a complementary measurement to the dileptons. 

3. Charge fluctuations 

Measuring the charge fluctuations or more precisely the charge fluctuations per unit degree 
of freedom of the system created in a heavy ion collision would tell us immediately if we have 
created a system of quarks and gluons [ |16| ] (see also [|r7|]). The point is that in the QGP phase, 
the unit of charge is 1/3 while in the hadronic phase, the unit of charge is 1. The net charge, 
of course does not depend on such subtleties, but the fluctuation in the net charge depends on 
the squares of the charges and hence are strongly dependent on which phase it originates from. 
However, as discussed in the previous section, measuring the charge fluctuation itself is plagued 
by systematic uncertainties such as volume fluctuations, which can be avoided if one considers 
ratio fluctuations. The task is then to find a suitable ratio whose fluctuation is easy to measure 
and simply related to the net charge fluctuation. 

The obvious candidate is the ratio F = Q/N c h where Q = N+ — N- is the net charge and 
A^ c h = N + +N- is the total charge multiplicity. Instead of using F, however, it is simpler to 
consider the the charge ratio R = N+/N-. If (A^h) 3> (Q), i.e. \F\ <C 1 the fluctuations of F are 
related to those of R by 



Furthermore in the limit of (Q) << (N c h) the fluctuations are dominated by the small term and 
we find [|T§ 



Usually, the number of charged particles is directly related to the entropy generated in these 
collision. Thus the observable 



provides a measure of the charge fluctuations per unit entropy. 

In order to see how this observable differs between a hadronic system and a QGP let us 
compare the value for D in a pion gas and in a simple model of free quarks and gluons. 

In a pion gas, the fundamental degrees of freedom are of course pions. Hence, A^h = N % + + 
N K - and using thermal distributions and disregarding correlations the charge fluctuations are 
given by 




(7) 




(8) 




(9) 




(10) 



where w n = (8N%) / (N n ) is slightly bigger than 1 [|5]|J. Hence for a pion gas, 



D 



71— gas 
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(11) 
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In the presence of resonances, this value gets reduced by about 30 % due to the correlations 
introduced by the resonances, as discussed in section (|2.3[). 

For a thermal system of free quarks and gluons we have in the absence correlation 

<5£ 2 ) = Qlw u (N u ) + Q 2 d w d (N d ) (12) 

where Q q is the charges of the quarks and N q denotes the number of quarks and anti-quarks. 
The constant w q = (SA^) / (A^) is slightly smaller than 1 due to the fermionic nature of quarks. 

Relating the final charged particle multiplicity A^h to the number of primordial quarks and 
gluons is not as simple. Using entropy conservation one finds [ fT6| ] 



N ch ~-(N g + 1.2N u +1.2N d ) (13) 
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leading to 



D QGP ~ 3/4 (14) 

Actually the charge fluctuations ((8<2) 2 ) have been evaluated in lattice QCD along with the 
entropy density [|TT]]. Using these results one finds 

DLattice-QCD ^ 1 — 1-5 (15) 

where the uncertainty results from the uncertainty of relating the entropy to the number of 



charged particles in the final state. Actually the most recent lattice result Q12| ] for the charge 



fluctuations, which was obtained in the quenched approximation, is somewhat lower then the 



result of [jry]. 

But even using the larger value of D = 1.5 for the Quark Gluon Plasma, there is still a factor 
of 2 difference between a hadronic gas and the QGP, which should be measurable in experiment. 

4. Observational issues and model calculations 

The key question of course is, how can these reduced fluctuation be observed in the final 
state which consists of hadrons. Should one not expect that the fluctuations will be those of the 
hadron gas? The reason, why it should be possible to see the charge fluctuations of the initial 
QGP has to do with the fact that charge is a conserved quantity. Imagine one measures in each 
event the net charge in a given rapidity interval Ay such that 

l«Ay« Ay mflX (16) 

where Ay max is the width of the total charge distribution. If, as it is expected, strong longitudinal 
flow develops already in the QGP-phase, the number of charges inside the rapidity window Ay 
for a given event is essentially frozen in. And if Ay ^ 1 neither hadronization nor the subsequent 
collisions in the hadronic phase will be very effective to transport charges in and out of this 
rapidity window. Thus, the E-by-E charge-fluctuations measured at the end reflect those of the 



initial state, when the longitudinal flow is developed. Ref. [ |T8| ] arrives at the same conclusion on 
the basis of a Fokker-Planck type equation describing the relaxation of the charge fluctuation 
in a thermal environment. In this approach the relaxation time for the charge fluctuations is 
determined by the average rapidity shift of the charges in a collision and by the number of 
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collisions in the hadronic phase. These values can be easily extracted from a transport model 
and we find that for Ay > 3 the signal should survive the hadronic phase. The smearing of the 
signal due to hadronization is more difficult to estimate, however. But it is rather unlikely that 
hadronization, which is a soft process, should give rise to large rapidity shifts. 

Thus, in the limit of very large Ay our proposed signal should be robust. However, one should 
not make Ay too large, as eventually charge conservation becomes relevant. Actually in the limit 
that Ay ~ Ay max charge conservation requires that the fluctuations vanish. For moderate values 
of Ay, one can very successfully correct for the effect of charge conservation, and this needs 
to be done for a meaningful interpretation of any experimental data. This correction factor is 
given by [jljj 

<> 7 > 

{■Nek) total 

Also in case of a finite net charge, which may be due to non negligible baryon stopping, some 
additional, though small, correction needs to be applied [|T9|] . 



Thus the actual observable turns out to be 
D 

D = (19) 

C C 

In Fig. @ the importance of these corrections, in particular the charge conservation correction 
C y is demonstrated. There we compare the uncorrected observable D with the corrected observ- 
able D as a function of the width of the rapidity window based on a URQMD [ |20| ] simulation. 
The effect of the charge conservation is clearly visible. With increasing rapidity window the 
charge fluctuation are suppressed. Once the charge-conservation corrections are applied, the 
value for D remains constant at the value for a hadron gas of D ~ 3. This is to be expected 
for the URQMD model, which is of hadronic nature and does not have quark- gluon degrees of 
freedom. 

Also for small Ay the value is D ~ 4 before it drops to D ~ 3 for Ay > 1.5. This effect, which 
was already predicted in [Q], is simply due to the fact that the correlation introduced by the 
resonances gets lost if the acceptance window becomes too small. 

In Fig. (g), we also show the predictions of the URQMD model as a function of impact 
parameter for SPS as well as RHIC energies JT^]. For both energies hardly any centrality 
dependence is visible and the results agree within errors with the prediction for the hadron gas. 
Thus a measurement of D ~ 1 would clearly indicate the existence of a QGP in the system 
created in these collisions. 



5. Some additional remarks 

A similar argument as for the charge fluctuations also holds for the fluctuations of the baryon 
number, since in the QGP the quarks carry fractional baryon number [JT7|]. The observation of 
reduced baryon number fluctuations, however, is much more difficult since one has to measure 
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Figure 2. Charge fluctuations as a function of the size of the rapidity window(left) and central- 
ity (right). 



neutrons on an event by event basis. Also, so far, no obvious intensive observable has been 
proposed which would be sensitive to baryon number fluctuations. 

Certainly, a measurement of D in proton-proton collisions will be necessary in order to rule 
out that possible small fluctuations are simply due to the nature of the proton wavefunction. 

In order to be able to see the reduced charge fluctuations, they have to be imprinted into the 
hadronic phase. How can this be possible? As we have shown above, the charge fluctuations of 
a hadronic gas are three times larger than that of a QGP. The only way to reduce the fluctuations 
in a hadron gas is to introduce additional correlations (see eq. @). As we have argued above, 
neutral hadronic resonances are the obvious candidates, which would lead us to predict an 
enhancement of neutral resonances. However, if, in addition, we require that on average iso- 
spin should be conserved, increasing the number of neutral rho mesons would also imply an 
increase of the charged rho mesons and thus large charge fluctuations. This leaves us with the 
iso-scalar resonances, such as eta, omega, phi etc. We thus are lead to predict an enhancement 
of eta, omega and phi production, which can be readily observed via electromagnetic decay 
channels. An actual calculation based on this conjecture gives enhancement factors of the order 
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of 5. 



6. Conclusions 

We have discussed event-by-event fluctuations in heavy ion collisions. These fluctuations 
may provide useful information about the properties of the matter created in these collisions, as 
long as the 'trivial' volume fluctuations, inherent to heavy ion collisions, can be removed. We 
have argued that the fluctuations of particle ratios is not affected by volume fluctuations. 

In particular the fluctuations of the ratio of positively over negatively charged particles mea- 
sures the charge fluctuation per degree of freedom. Due to the fractional charge of the quarks, 
these are smaller in a QGP than in a hadronic system. 

A measurement of our observable D ~ 1 would provide strong evidence for the existence of 
a QGP in these collisions. A measurement of D ~ 3 on the other hand does not rule out the 
creation of a QGP. There are a number of caveats (see flT^]), which could destroy the signal, 
such as unexpected large rapidity shifts during hadronization. 
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